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THEORETICAL STUDIES OF SUPERSONIC TWO-DIMENSIONAL AND 
AXISYMMETRIC NONEQUILIBRIUM FLOW, INCLUDING 
CALCULATIONS OF FLOW THROUGH A NOZZLE 1 
By James J. Der 

SUMMARY 


Chemical and vibrational nonequilibrium phenomena in steady two-dimensional 
and axisymmetric inviscid flow fields are studied by the analysis of flow past 
curved boundaries. The study consists of three parts: (l) formulation of the 

governing equations, (2) study of the general features of nonequilibrium flow by 
consideration of some mathematical properties of the nonlinear flow equations, 
and ( 3 ) analysis of axisymmetric flow of dissociated air through a nozzle by 
means of numerical computations. 

In the formulation of the governing equations, possible alternative choices 
of the state variables to be used and their relative advantages are considered. 
The entropy is not found to be a convenient variable in nonequilibrium flow and 
is therefore not used. The system of nonequilibrium flow equations is two orders 
higher than the corresponding classical (equilibrium or frozen) system. 

General features of the nonequilibrium flow are studied by examination of 
the mathematical properties of the nonlinear flow equations. In addition to the 
two differential compatibility equations along the Mach lines, as exist in clas- 
sical supersonic flow, two more differential compatibility equations are neces- 
sary along the streamlines. The domain of determinacy of the nonequilibrium flow 
is the same as for classical flow, namely, the area bounded by the outermost Mach 
lines from a data curve. If a data curve is crossed more than once by any char- 
acteristic, the amount of data that can be prescribed is restricted. The 
implications of such restrictions in nozzle problems are discussed. 

Finally, an axisymmetric nonequilibrium flow of air through a supersonic 
nozzle is analyzed by means of numerical computations. The coupling effects 
between nonuniformity (across a nozzle section) and nonequilibrium, not taken 
into account in the quasi -one -dimensional flow approximation, are revealed by 
the present analysis. 

iThis report is based on a dissertation submitted to the Department of — 
Aeronautics and Astronautics, Stanford University, in partial fulfillment of the 
requirements for the Ph.D. degree. The work arose out of the author* s participa- 
tion in a graduate research seminar supported at Stanford University by a grant 
from the National Science Foundation. 


CHAPTER I 


INTRODUCTION 


The present study is concerned with chemical and vibrational nonequilibrium 2 
phenomena in two-dimensional and axisymmetric flow. 

Studies of nonequilibrium flow have been concerned mainly with one- 
dimensional (refs. 1 through 3) and quasi -one -dimensional flow (refs. 4 
through 11) and with two-dimensional and axisymmetric flow past bodies of simple 
shape (refs. 12 through 17 ). 3 A general formulation of the equations for non- 
equilibrium flow has been presented by Kirkwood and Wood (ref. 18 ) and by others \ 
(e.g., refs. 19 and 20). A method for the numerical computation of the quasi- 
one -dimensional flow of a general gas mixture has been formulated recently by 
Emanuel and Vincent! (ref. 8). 

In the present work, the more general features of two-dimensional and 
axisymmetric nonequilibrium flow are studied by consideration of flow past curved 
boundaries. For simplicity the study is limited to the cases of supersonic flow 
with a single nonequilibrium process. Part of this objective was achieved in 
reference 21, which, however, was concerned only with linearized flow past a 
single wall. 

In chapter II, the inviscid flow equations are given. Alternative forms of 
the governing equations and their possible advantages are discussed. The 
entropy, a convenient quantity in classical (frozen or equilibrium) gas dynamics, 
has no particular utility in the present analysis as a none qui librium variable. 

The constancy of a nonequilibrium parameter (other flow quantities varying) 
signifies the freezing of the corresponding nonequilibrium process. On the other 
hand, the entropy is constant along a streamline when the flow is either in an 
equilibrium or a frozen condition. The entropy is, therefore, not used in the 
formulation of the flow equations in the present work. The system of differen- 
tial equations for nonequilibrium flows is seen to be two orders higher than the 
corresponding system for classical gas flows. 

In chapter III, the general features of steady two-dimensional and 
axisymmetric nonequilibrium flow are studied by examination of some mathematical 
properties of the nonlinear flow equations. As in the case of classical gas 
dynamics, there are two differential compatibility relations along characteris- 
tics having the direction of the frozen Mach lines. In addition, however, there 
are two more differential compatibility relations necessary along the character- 
istic in the streamline direction. The domain of determinacy of the nonequilib- 
rium flow is the same as that of the frozen flow, namely, the region bounded by 
the outermost frozen Mach lines from a data curve. 


2 When a nonequilibrium variable (such as a degree of dissociation) is in 
equilibrium with the other local thermodynamic state variables, the flow is said 
to be in equilibrium with respect to the corresponding nonequilibrium process; 
otherwise the flow is said to be in nonequilibrium condition. When a nonequilib- 
rium variable becomes fixed, the flow is said to be frozen with respect to the 
corresponding nonequilibrium process. 

3 The references cited here constitute a representative rather than an 
exhaustive list of the numerous published works in this field. 
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For any system of hyperbolic differential equations, if a data curve is 
crossed by any characteristic line only once, full data can be posed on the data 
curve. If a data curve is crossed by any characteristic twice, then only partial 
initial values can be prescribed. For the system of equations used herein, the 
prescription of initial values of the pressure and stream angle will be 
restricted if the data curve is crossed twice by one of the frozen Mach lines 
(which is true also for the classical flow) . On the other hand, prescription of 
the initial values of the velocity and the degree of dissociation will be 
restricted if the data curve is crossed twice by the streamlines. Nonequilibrium 
flow has more compatibility relations than the classical flow and hence has more 
restrictions on the prescription of initial values. These additional restric- 
tions prevent us from prescribing a completely uniform exit flow in the inverse 
nozzle problem (i.e., computation of the contour of a nozzle for given initial, 
center line, and exit conditions). 

Finally, a specific nonlinear, nonequilibrium flow problem is studied in 
chapter IV, where the dissociative nonequilibrium flow of air through a given 
nozzle is analyzed. The flow of real air is approximated by a flow of a mixture 
of dissociating oxygen and nondissociating nitrogen. The molecular vibration is 
assumed to be in equilibrium. The assumption of only a single nonequilibrium 
process simplifies the analysis considerably. For the present nozzle problem 
(temperature less than 8000° K, pressure less than 300 atmospheres), however, 
such approximation is valid since nitrogen dissociation is actually small. Four 
cases for the same reservoir conditions were computed numerically by means of an 
IBM 7090 digital computer. These include one equilibrium flow and three cases of 
nonequilibrium flow with varying degrees of nonequilibrium as dictated by the 
magnitude of the effective characteristic reaction length. 

It is found that quasi -one -dimensional flow computations give qualitatively 
correct results. On the other hand, the coupling effects between nonunif ormity 
(at a nozzle cross section) and nonequilibrium, not taken into account in quasi- 
one -dimensional flow, are revealed in the present analysis. At any cross sec- 
tion, the flow near the nozzle wall has traveled a longer distance than the flow 
near the nozzle center. Furthermore, the flow near the center line freezes 
sooner. Hence, near the center line the flow has higher values of the degree of 
dissociation. The nonuniformity at the nozzle exit worsens as the degree of non- 
equilibrium increases. The exit velocity changes only minutely and remains 
nearly uniform as the flow deviates from the equilibrium condition. Other flow 
variables, including pressure, temperature, and Mach number, change significantly 
both in magnitude and uniformity at the exit as the flow tends toward the frozen 
condition. In general, therefore, one can expect little effect on velocity due 
to nonequilibrium, but significant effects on Mach number and other variables. 
Also, since one cannot impose complete uniformity of all flow variables at the 
nozzle exit, prescribing uniform exit pressure in designing a nozzle would lead 
to more uniform flow over -all than prescribing uniform exit velocity. 

The thermodynamic and chemical kinetic properties of the simplified air 
model used in the numerical analysis are listed in appendix A. 

The author is indebted to Professor Walter G. Vincent! of Stanford 
University for his guidance throughout the course of the research and to 
Dr. Max. A. Heaslet of the National Aeronautics and Space Administration, Ames 
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Research Center, for his timely encouragement and assistance. Professor 
M. D. Van Dyke provided valuable criticism and encouragement . Finally, thanks 
are due to Professor I -Dee Chang for reading the manuscript and furnishing 
important criticism and assistance. 
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CHAPTER II 


FLOW EQUATIONS 


The governing equations of flow of a reacting gas in the absence of 
transport phenomena and body forces are the familiar conservation equations plus 
the equation of state and the rate equations. The rate equations, which are not 
used in classical (equilibrium or frozen) gas dynamics, are needed here because 
the more general state equation now contains additional nonequilibrium param- 
eters, for example, the various molecular vibrational energies and the degree of 
dissociation and ionization of the various constituents. The equations for flow 
with n nonequilibrium parameters are, therefore (see e.g., ref. 22) 


Mass 

Momentum 

Energy 

State 

Rate 


?£ + p div w = 0 
Dt 


(II -la) 

p .— + grad p = 0 
Dt 


(II -lb) 

p “i _ 5E = 0 
Dt Dt 


(Il-lc) 

h = h(p, p, q_ ± , ... 

' %) 

(II -Id) 

Dqt 

_ = UiCp, P, q x , ... 

• ) On) 

i = 1, ... , n 

(Il-le) 


where D/Dt denotes the Eulerian derivative; p, w, p, h are the density, 
velocity, pressure, and enthalpy of the medium; q^ denotes the vibrational 
energy and/or degree of dissociation and ionization; and gives the net rate 

of increase of q^ as a function of the local thermo chemical state. 


When the flowing gas is in equilibrium, the values of q^_ are functions 
solely of the local pressure and density, hence the enthalpy becomes a function 
solely of p and p, that is 


h e = h[p, p, q. ie (p, p)] = h e (p, p) 

When the gas is in a frozen condition, each q^ has a certain fixed value, and 
h is again a function solely of p and p, that is, 


hf = h[p, p, q.j_ = constant] = hf (p, p) 

Note that a perfect gas has qj_ = qp (p, p) = constant; hence a perfect gas is, 
by definition, always in the frozen and equilibrium conditions simultaneously. 
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It is of interest to note the minimum number of first-order differential 
equations required to determine a flow field. We consider, for simplicity, 
steady two-dimensional and axisymmetric flow with only one mode of nonequilib- 
rium. Since the state equation is not in differential form, it can be used to 
eliminate one dependent variable from the differential equations. We shall 
choose to eliminate p. This can be done by differentiating equation (il-ld) 
and using equations (il-lc) and (il-le) together with the identity 



(II -2) 


for the frozen speed of sound (see appendix A), thus finding an expression for 
Dp/Dt in terms of Dp/Dt and Dq^/Dt. By substituting this expression into equa- 
tion (il-la), we obtain a modified mass equation. Furthermore, the energy equa- 
tion can be modified by combining it with the momentum equation, yielding a form 
more convenient for our purpose. A complete set of differential equations 
governing nonequilibrium flow is thus obtained as 


ii 

“ 2 2E + p div w - ^ 


JD_ 

dt 


n V 

V h<1 i 

• ) -^“1 = 0 
1=1 

(II -3a) 

— + grad p = 0 
Dt 

(H-31) 

L + = 0 

2 J <3t 

(II -3c) 

DcLj. 

1 - ( j -: - 0 

Dt 1 

(II -3d) 


The coefficients are related to the dependent variables, p, w, h, q^. 

For steady flow, equation (ll-3c) can be integrated immediately, yielding an 
algebraic relation stating that the sum h + w^/2 is constant along each stream- 
line. We further restrict our study to two-dimensional axisymmetric flow, with 
only one nonequilibrium parameter. Equations (II-3) then become, in intrinsic 
coordinates. 


3a 


Q S vw ) = 0 

hp J 

(II -4a) 

pw = o 

3s 3s 

(II -4b) 


r 



(II -4c) 


pw s + = o 

bs bn 


- - = o (ii -4a) 

5s w 


where Mp = w/a^; s and n are the streamwise and normal coordinates; 6 is the 
streamline angle; v is 1 for planar flow and 2 for axisymmetric flow; and r is 
the radial distance from the center line of axial symmetry. 


The minimum number of first-order differential equations governing a steady, 
nonequilibrium, two-dimensional or axisymmetric flow is thus seen to be four. 


In the classical limits of either equilibrium or frozen flow, two more 
differential equations can be eliminated. This can be done by introducing the 
specific entropy S, which is related to by the relation 


T 


DS 

Dt 


n 



i=i 


(II-5) 


where the z^ are related to the chemical potential and the vibrational energies 
(see, e.g., ref. 18) . (The specific entropy is denoted here by S to avoid con- 
fusion with the streamwise coordinate s.) The rate of change of entropy DS/Dt 
is zero when the flow is frozen (up = 0 ) or in equilibrium (]^ zpwp - 0 ). Fur- 
thermore, since the state of the medium in either limit is determined by only two 
thermodynamic variables, we can take one of these to be the entropy, which is, 
from equation (II-5), constant along streamlines. The two relations, h = h(p, S) 
and S = constant along a streamline, can then be used to replace two of the dif- 
ferential equations, equations (ll-4b) and (ll-4d) . In the frozen and equilib- 
rium limits, therefore, the order of the system of flow equations is 2 . This is 
true regardless of whether the flow is rotational or irrotational . 

The governing differential flow equations can be written in various 
alternative forms. For special purposes a particular form of the equations may 
have advantages over the others. The differences among various forms of the 

equations arise mostly from the choice of the state equation. Choice of equa- 
tion (il-ld) is natural because it relates h, which arises naturally in the 

energy equation for fluid flow, to p and p, which appear in the mass and 

momentum equations. Equations (il-l) are convenient for analytical purposes. 

A perhaps more familiar way of expressing the state relationship is to 
relate h, p, and p by means of the two equations 


h = h(p, T, q. i ) 

(II -6a) 

P = p(p, T, q. ± ) 

( II -6b ) 
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•with the temperature T as a parametric quantity. These two equations are 
called, respectively, the caloric and thermal equations of state. For numerical 
computations, equations (II-6) are often convenient to use because the tempera- 
ture usually appears as a natural variable in the thermodynamic and kinetic 
relations for real gases, such as air (see appendix A). 

The entropy, useful in classical gas dynamics, especially when the flow is 
homoentropic, is not nearly so convenient in nonequilibrium flow since here 
entropy is not constant along streamlines. On the other hand, the constancy of 
the q^ (with other flow quantities varying) signifies that the flow is in the 
frozen condition. Furthermore, the derivative of entropy in equation (II-5) con- 
tains Zj_ and Gdq, and in some conditions both of these become indeterminate. The 
use of entropy, therefore, may affect the stability of the numerical integration 
process. If one needs to examine the variation of entropy, it is preferable to 
compute it from equation (II-5) after the flow has been found in terms of other 
variables . 
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CHAPTER III 


SOME MATHEmTICAL PROPERTIES OF THE NOl^QUILIBRIUM-FLOW EQUATIONS 


We no w consider the application of the flow equations (eqs. II-4) developed 
in the last chapter to problems of fluid flow in nonequilibrium condition. These 
equations are quasi -linear 4 with complicated coefficients. Practical solutions 
to such equations can be obtained, except in special cases, only by numerical 
methods. We can study some of the general characteristics of nonequilibrium 
flow, however, by examining a few mathematical properties of the nonequilibrium 
flow equations and comparing them with those of the classical flow equations. 

For simplicity we restrict our study to cases of only one nonequilibrium 
process since the extension to multiprocess cases introduces no additional dif- 
ficulties. Again, we study only steady two-dimensional or axisymmetric flow. 
Supersonic flow will be dealt with mostly, although some of the discussion 
applies to more general speed regimes. After some general discussion of the 
mathematical properties, we apply the results to nozzle problems with regard to 
the posing of initial and boundary values. Strong discontinuities such as shock 
waves and slip lines will not be considered. 

Some mathematical questions of special Interest are the following: 

1. What is the minimum number of dependent variables that determine the 
flow field? 

2. What are the characteristics? 

3* What must be continuous, and what can be discontinuous, across the 
characteristics? 

k. What is the domain of determinacy? 

5* How can the initial values be prescribed? 

6. What can be said about the existence and uniqueness of the solution? 

We may try to study these questions in terms of either a single nth order 
differential equation or a system of n first-order equations. We choose the 
latter because, as pointed out by Courant (ref. 23 ): 

l. A higher order differential equation can be rewritten as a system of 
first-order equations, but the reverse is not true in general. 

2. It is clearer to speak of the fundamental” flow quantities, such as 
velocity and pressure, than of some relatively remote function of the flow 
variables, such as the stream function. 


4 A quasi-linear differential equation Is linear in its highest derivatives. 
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3 . The gas-dynamic equations can he -written readily as a system of first- 
order equations. 

The minimum number of dependent variables that determine the flow field is, 
of course, the same as the minimum number of differential equations. As shown 
in chapter II, the minimum number of differential equations in the present case 
is four. The dependent variables selected are p, 0, w, and q, but any 
appropriate function of these variables could also be used. 

The characteristics of a system of differential equations can be defined in 
various ways. The definition used by von Mises (ref. 24) is as follows: 

"Curves along which analytically different solutions of a 
differential equation or a system of such equations can be patched 
together are called characteristics of the equation or system." 

Here patching means the solutions on both sides of the curve are equal. Charac- 
teristics are thus curves across which the solutions themselves must be contin- 
uous but the normal derivatives of any order can be discontinuous. This type of 
discontinuity is called a weak discontinuity. 

An equivalent definition of the characteristics is that along a characteris- 
tic certain partial derivatives of the equations group together in such a manner 
that they form a differential equation containing only ordinary derivatives in 
the direction of the characteristic. The characteristics and compatibility 
relations for the nonequilibrium -flow equations can readily be obtained by using 
this last definition. The direction along which an ordinary differential equa- 
tion is formed Is then the direction of a characteristic. The ordinary differen- 
tial equation itself provides the compatibility relation corresponding to this 
characteristic. The compatibility relations are, in fact, the original differen- 
tial equations written in characteristic form. The order of the system remains 
unchanged. Thus, in the present case we look for four compatibility relations. 

Equations (ll-4b) and (ll-4d) are already in characteristic form along the 
streamlines. After elimination of the derivative of w by means of equa- 
tion (ll-4b), the partial derivatives in equations (ll-4a) and (ll-4c) can be 
grouped in such a manner that these two equations form a set of two ordinary 
differential equations in terms of p and 9 along the two frozen Mach lines. 
Briefly, these equations can be obtained as follows: If we eliminate dw/ds by 

the use of equation (ll-4b), equation (ll-4a) becomes 

(M-p 2 - l) + pw 2 — + ( v ~ 1 pw 2 sin 6 wto') = 0 

Ss Sn \ r hp / 

The total differentials of p and 6 in the Z direction (defined by the slope 
dn/ds of Z) can be written as 

dp = ds + ^2 dn 
Ss 3n 
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and 


ae = — ds + dn 
bs on 


Solving these equations for bj>/bs and bd/bn in terms of the directional deriv- 
atives dp/ds and d0/dn, substituting into the foregoing equation, and noting 
that p\r 2 bo/bs = -Sp/5n (from eq. (ll-4c)), we have 


Pf" 


dp 

ds 


P W2 M + 
dn 


ds • q 2 . 

Pf 

dn 


dn 

ds. 


ap 

bn 



sin Q 



= 0 


If we choose dn/ds = ±l/Pf, then the term containing the partial derivative 
bp/bn vanishes, and we thus obtain two ordinary differential equations along 
Z x and Z^, that is, along 


dn 

ds 


1 

p f 


or 


dZ x _ Mf 

Pf 


ds 


and 


dn 

ds 


_1_ 

Pf 


or 


dZ 


ds 


2 _ 


Mf 

P f 


In terms of Zi, 1%, therefore, we have 


Pf* dp ± pw 2 d0 + ( v - ~ 1 pw 2 sin 0 - — wto ) dZi 2 = 0 
1 v Mfr Mfhp J 9 


\ 


(ill-la/b) 


and 


along 


— = tan(0 ± \i) 
dx 


pw dw + dp = 0 
w dq_ - co ds = 0 




along 


^ = tan 0 
dx 


(ill-lc) 
(ill -Id) 
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where = JMp 2 - 1; l x and are the distances along the left-running and 

right-running Mach lines, respectively (fig. III-l); r is the radial distance 
from the line of symmetry, which is also the x-axis; and \i is the Mach angle, 
related to Mf "by the relation |u = sin _1 (l/Mf ) . 


r 



Figure III-l.- Characteristic lines for equations (ill-l) . 

There are other, alternative forms of the compatibility relations. They can 
he obtained by the use of different variables. In general, the number of com- 
patibility relations needed is the same as the number of original differential 
equations. Normally, there are two along the slant characteristics, or Mach 
lines, and the remainder along the streamlines. In selecting a particular set of 
compatibility relations, or variables, one should take into account the manner in 
which the problem is prescribed. 

We can now deduce from equations (ill-l) the answers to the questions posed 
earlier. First of all, p, w, 0, q constitute the solution to the system of 
equations, and hence these flow quantities must be continuous everywhere in the 
domain of weak discontinuities (i.e., in a flow field free of shock waves and 
slip lines). Equation (ill-la) is the characteristic equation for p and 0 
along the left-running Mach line, whereas equation (ill -lb) is the characteristic 
equation for p and 6 along the right-running Mach line. Across any Mach line, 
therefore, the normal derivatives of p and 0 can be discontinuous. Equa- 
tion (ill-lc), on the other hand, is the characteristic equation for p and w 
along the streamlines. It follows that the normal derivatives of p and w can 
be discontinuous across a streamline. Similarly, due to equation (ill-id), the 
normal derivative of q can be discontinuous across a streamline. 
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Note that, alternatively, p instead of w could have "been eliminated in 
forming equations (lll-la,b) . The variable w, therefore, can have a weak dis- 
continuity (normal derivative discontinuous) across the Mach lines also. The 
quantity 0, however, cannot have a weak discontinuity across a streamline. This 
can he deduced by the use of equation (ll-4a) as follows: Since all of the 

coefficients and, in general, 3p/3s and 3w/3s are continuous across the stream- 
line, It follows that 30/3n must be continuous across a streamline. Similarly, 
that the quantity q cannot have a weak discontinuity across a Mach line can be 
deduced from the fact that 3q/3s is always continuous (since w is always 
continuous) and dq/ 37 i ,2 are in general continuous (since q is in general 
continuous ) . Summar iz ing : 

p, w, 0 can have weak discontinuities across a Mach line, while 
p, w, q can have weak discontinuities across a streamline. Weak 
discontinuities cannot occur in q across Mach lines nor in 0 
across streamlines. The actual occurrence of any such disconti- 
nuity, of course, is dependent on the initial conditions. 

From the point of view of analysis, the characteristic equations have two 
properties of significance. One is that they are ordinary differential equations 
along the characteristic lines, and hence enable us to use a relatively well 
developed numerical integration scheme. (This scheme is discussed in chap- 
ter IV.) The other is that the variation of the flow quantities along the char- 
acteristics cannot be arbitrary but is restricted by the relationships fixed by 
the characteristic equations . These restrictions limit the freedom with which 
initial values of p, w, 0, q can be prescribed. The nonequilibrium flow 
requires a larger number of characteristic equations than the classical flow and 
hence has more restrictions . We can more readily visualize the problem of posing 
the initial conditions after the domain of determinacy is discussed. 


The theory of characteristics has been studied extensively for linear 
equations. The results of these studies can be extended by successive approxi- 
mations to quasi-linear equations. Here, we are specifically interested in the 
domain of determinacy. Consider figure III -2 where initial data are prescribed 



'Characteristics 


Figure III-2.- Domain of determinacy for equations (ill -la, b) . 
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on the data curve AB. By data curve we mean a segment of a curve along which 
the initial values of the dependent variables of a system of differential equa- 
tions are prescribed. If the initial values of all the dependent variables of 
the system of equations are prescribed, this segment is said to have full data, 
otherwise the data curve is said to have partial data. A data curve may or may 
not constitute the data curve of the entire problem (flow field in the present 
case); and the system of equations may be either the complete system, that is, 
all of equations (ill-l), or a subsystem of it, for example, equations (ill -la, b). 

If the subsystem of equations (ill -la, b) were linear, the domain of 
determinacy, with full data of this subsystem prescribed on AB, is the area 
ACBD bounded by the outermost slant characteristics (frozen Mach lines) AC, BC, 
AD, and BD. 

For our equation, which is quasi-linear instead of strictly linear, we can 
apply the iterative procedure of Schauder (ref. 23, P* 476) to compute the flow 
quantities to any desired accuracy in the small region A’B^*, where A ! C T and 
B ! C ! are characteristics. By continuing the process over the full length of 
AB, we can determine the flow quantities in the region AA"B ,T B in the near 
neighborhood of AB. We thus obtain a new data line A tf B". By continuing the 
process toward C and also toward D, we can determine the entire region ACBD 
bounded by the outermost characteristics. For the nonlinear equations, of 
course, the location of the characteristics is not known in advance, but must be 
computed successively. According to reference 23 (p. 466), the above Iterative 
solution converges and is unique, and the solution can be extended to a finite 
domain containing weak discontinuities. Hence, the solution to the quasi-linear 
equations exists and is unique in domain ACBD. 

In our system of equations, we have not only slant characteristics, that is, 
the Mach lines, but also the streamwise characteristics defined by the stream- 
lines. If equations (lll-lc,d) were linear, the domain of determinacy of each of 
these subsystems (eq. III-lc, eq. Ill-Id) would be an infinite strip bounded on 
the sides by the outermost streamlines passing through the initial data line AB 
(fig. Ill -3)- 



Figure III-3-- Domain of determinacy for equations ( III-lc, d) . 
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All of the subsystems (eqs. Ill -la, b; III-lc; III -Id) are, of course, 
coupled. For example, the coefficients for equations (lll-la,b) are functions 
not only of p and 9 hut also of w and q. Equation ( III-lc) likewise has a 
coefficient that is a function of p, v, and q, and, in addition, the quantity 
p appearing in the differential is to he found from the subsystem equa- 
tions (ill -la, h). The domain of determinacy of the complete system of equa- 
tions (ill-l) is, therefore, the common domain of determinacy of all the sub- 
systems. Thus, the domain of determinacy of the complete system is the same as 
that for subsystem equations (lll-la,b), that is, the region bounded by the 
outermost slant characteristics. 

In posing the initial conditions we must take care that no region is 
overdetermined. In specifying the Initial values of p and 0, for example, we 
must inquire whether the data curve is crossed more than once by any Mach line 
(since the Mach lines are the characteristics for the subsystem equa- 
tions (lll-la,b)). To see this we refer to figure III - 4 . The segments AB, BB 
and BB^ are data curves. If the initial values of p and 9 (as well as w 
and q) are prescribed on AB, then, since BBi lies Inside the domain of deter 
minacy of AB, the values of p and 0 on BB]_ are determined. Hence, neither 
nor 9 can be prescribed on BB X . In other words, if a curve (e.g., ABBi) is 



Figure III-4.- Data curves with overdeterminacy . 





crossed twice by both the left-running and right-running Mach lines, then one 
segment of it (BB X ) lies inside the domain of determinacy (for AB) of the sub- 
system equations (ill-la^b). Consequently, the full data of this subsystem 
(values of p and B) can be prescribed only on one segment (AB) . Posing partial 
data on both segments (AB and BB X ) is possible, however. For example, if p 
only is prescribed on AB and 0 only is prescribed on BB X (and also on 
BiB 3 ), then, p for BB X and Q for AB are determined by equations (ill -la, b). 

Another situation that can arise is illustrated as follows: The curve 

A 2 BB 2 is crossed twice by only one of the two families of Mach lines. The data 
prescription on such a curve is less restrictive than that on a curve of the 
type ABB X . If full data for the subsystem equations (lll-la,b) (both p and 6) 
are prescribed on A 2 B, then the values of p and Q on BB 2 are related by equa- 
tion (ill-lb). Consequently, either p or 0, but not both, can be prescribed 
on BB 2 . 


In posing data for w and q, we need to note whether the data curve for 
w or q is crossed more than once by a streamline. In figure III -4, AiBBi is 
crossed twice by the streamlines. If the initial value of w is prescribed on 
Aj_B then the value of w on BB X is related to that on A X B by equa- 
tion (ill-lc) . The initial value of w (and similarly q) can be prescribed 
on either A X B or BB X but not on both. 

The foregoing are the restrictions imposed by the characteristic relation 
as mentioned earlier in this chapter. We can summarize these restrictions as 
follows (referring to fig. III-4) : 

I. If a data line is crossed twice by 

A. both families of the Mach lines (say first segment AB, then 

segment BB X ) , then the initial value of p and B can be posed 
in one of two ways : 


1. 

Prescribe 

p 

and. 0 

on AB 

or BBi 

only. 

2. 

Prescribe 

p 

on AB 

and 0 

on BBjl 

or vice versa 


B. only one family of Mach lines, then both p and B can be 

prescribed on one of the segments and either p or B on the 
other segment. 

II. If a data line is crossed by the streamlines twice (say first A X B, 
then BB X ), then w can be prescribed on one of the two segments 
only. The same is true for q. 

We illustrate these points by the following two examples. 
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Direct Nozzle Problem 


In the direct nozzle problem, for which calculations are made in the next 
chapter, initial data can be prescribed on AB (fig. Ill -5) • For our system. 



Figure III-5-- Direct nozzle problem. 

the data that can be prescribed on this line are p, w, 0, q.. The values of 0 
are also prescribed on AD where 0 = 0 and BC where 0 = 0(x) . Curve BAD is 
crossed twice by only the right -running Mach lines and therefore p and 0 can be 
prescribed on AB and 0 alone can be prescribed on AD. Similarly, 0 alone 
can be prescribed on BC . Curve ABCD is crossed by the streamlines twice; 
hence w and q can be prescribed either on AB or CD, but not both. 

The posing of initial values in the direct nozzle problem is thus straight- 
forward even for nonequilibrium flow. The inverse nozzle problem, on the other 
hand, is more deceptive. 


Inverse Nozzle Problem 


In the inverse nozzle problem, usually the initial values are given in the 
throat region; certain desired exit and center-line conditions are also 
prescribed. The problem is to find the nozzle contour that gives these 
conditions . 
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Suppose ve prescribe the initial values of p, 0, w, q on AB (fig. Ill -6) 
which is a curve crossed only once by any characteristic. The value of 0 = 0 is 



Figure III-6.- Inverse nozzle problem. 


prescribed on AD by the definition of the nozzle center line. On the center 
line, therefore, ve can prescribe p only on A T D ! since both BAA f and D r DC 
are crossed twice by a family of the slant characteristics. On CD, which is a 
curve crossed also only once by any characteristic, the values of p and 0 can 
be prescribed; but w and q cannot be prescribed on CD because ABCD is 
crossed twice by the streamlines and the initial values of w and q are already 
prescribed on AB. Summarizing, we can prescribe all of the flow quantities 
p, 0, w, q on AB; 0 on AD; p on A T D T ; and p, 0 on CD. Alternatively, A f B and CD r 
can be used as data lines in place of BAA 1 and CDD T , respectively. Then all of 
the flow quantities p, 0, w, q can be prescribed on A T B, and p and 0 can be 
prescribed on CD ! . 

In nozzle design it is usually desired that the flow at the exit be uniform. 
If CD in figure III -6 is perpendicular to the nozzle center line, we can let CD 
represent the exit plane. To prescribe that the flow be parallel on the exit 
plane, we let 0 = 0 on CD. To have completely uniform flow we would further 
require that p, w, q be constant along CD. However, as previously concluded, 
we can prescribe only p and 0 on CD. We can therefore have only a parallel flow 
on the exit plane with p uniform, with no assurance that w and q will also be 
uniform. Alternatively, we could, of course, prescribe w and q as uniform on 
CD instead of having given values on AB. Unfortunately, however, prescribing 
w and q on CD may dictate unrealistic conditions at the entrance AB of the 
nozzle (the same is true, of course, for the direct nozzle problem). 
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For the design of a nozzle with the flow either frozen or in equilibrium, 
the above dilemma does not occur because in these cases the specification that 6 
be zero and p uniform on CD leads automatically to uniformity of all other 
flow quantities on CD. 

Finally, we recall that w can exchange roles with p in the formulation 
of the characteristic equations. We can therefore alternatively prescribe w 
instead of p to be uniform on CD (q being unspecified in both cases). The 
results of chapter IV will show, however, that w is nearly uniform anyway for 
a reasonably shaped nozzle. As a practical matter, therefore, the prescription 
of uniform p rather than w on the exit is to be preferred. 
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CHAPTER IV 


AXISYMMETRIC NONEQUILIBRIUM FLOW THROUGH A NOZZLE 


We now study the flow of air through a high-speed nozzle where the air is 
out of chemical equilibrium. Specifically, we study the physical phenomena 
associated with dissociation nonequilibrium when the flow is not uniform across 
each section of the nozzle. We can thus examine the interrelated effects of 
nonequilibrium and nonuniformity. 

A typical nozzle is sketched in figure IV -1. To have significant effects 
of nonequilibrium, the air in the reservoir should be substantially dissociated, 
which means that the reservoir temperature must be high. This is usually the 
condition in a shock tunnel or spark -heated tunnel. 



Figure IV-1.- Regions of flow in a typical high-speed nozzle. 

In the region near the reservoir the high temperature of the air gives it a 
small relaxation time: In view of its small velocity, the air thus has an 

extremely short relaxation length here. In this region, therefore, the flow is 
in many cases in equilibrium. As the air approaches the throat the relaxation 
time greatly increases due to the dropping of the temperature. Also, the speed 
rapidly increases from subsonic toward supersonic values. Often, therefore, the 
relaxation length becomes moderate in the vicinity of the throat, and, if the 
relaxation length becomes moderate before the dissociated atoms have recombined, 
nonequilibrium conditions exist. As the air continues downstream, the relaxation 
length becomes large, and a near-frozen condition may exist at the exit of the 
nozzle . 


Gas Model 


In the preceding work to keep the discussion general, we have not been 
explicit as to the particular type of nonequilibrium, nor have we specified the 
actual equation of state. Now, to obtain a particular solution by numerical 
methods, we must decide on the gas model. A gas model that is a good approxima- 
tion to real air for a moderate range of temperatures, and yet relatively simple, 
is a mixture of dissociating oxygen and nondissociating nitrogen. The molecular 
vibration is assumed to be at equilibrium, as are the molecular translation and 
rotation. The expressions for the thermodynamic and chemical -kinetic properties 
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of this simplified air are briefly derived in appendix A. In this derivation, 
the contribution of electronic excitations to the internal energy is neglected. 
This is equivalent to assuming the electronic partition functions to be constant, 
which is valid for the temperature range concerned (T < 8000° K) . 

The difference of the present air model from the Lighthill ideal dissociat- 
ing gas (refs. 25 and 26 ) is due mainly to the inclusion of the dissociatively 
inert nitrogen. The principal approximation made in the Lighthill gas model, 
that pjj (defined by eq. (A-5)) is constant, is rather good for the range of 
temperature where dissociation is important. The Lighthill gas model, however, 
is not a good approximation for air for the present purpose because it has too 
high a value of ap/a e , which is an important parameter in estimating the non- 
equilibrium effects. The use of too high a value of a^/a e would lead to an 
overestimate of the effects of nonequilibrium. 

The simplicity of the present air model, as compared with more complex ones, 
is due mainly to the fact that there is only one nonequilibrium process. For the 
higher temperature range where dissociation of nitrogen is important, an equally 
simple but not as good approximation for air would be a mixture of dissociating 
nitrogen and fully dissociated oxygen. 


Nozzle Model 


The dimensions of the nozzle used in the present computations are shown in 
figure IV -2. The nozzle has a conical portion after the throat, followed by a 
fourth-degree -polynomial curved portion that matches the slope at the end of the 
conical section and has zero slope and curvature at x = 10. 35 cm. The conical 
section is included because a method of computation for the transonic throat 


Notes 

1 Curved portion at nozzle described by 

radius = 0.21256 x - 0. 01 8783 x 2 + 0 0001 42 2 6 * 3 + 0 000031201 * 4 

2 Throat at Jr = 0 2352 cm t radius at throat=0. 05cm 

3 Expansion area ratio at exit=387.75 



Figure IV-2.- Nozzle used in flow computations. 
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region is not available. If the conical section is large enough compared with 
the throat region, however, it is reasonable to assume that the flow itself is 
conical at the end of such a section. We select a fourth -degree -polynomial 
contour so that we can specify that the conical and the curved portions of the 
nozzle have matching slopes at their junction and that the curved portion has 
zero slope and curvature a short distance downstream from the nozzle exit. The 
selection of the location of the exit is described in the next paragraph. Such 
a nozzle may be thought of as representing the core of inviscid flow in a high- 
speed wind tunnel when the flow behavior differs (in the inviscid core or in the 
viscous layer, or both) from the design condition of uniform exit flow. In 
reality, the presence of a thick boundary layer in a wind-tunnel nozzle, such as 
the one we are considering, makes it difficult to achieve an exactly uniform flow 
in the test section under any condition, equilibrium or otherwise. 

Since this nozzle is not ideally designed, that is, its wall is not so 
contoured that the exit flow is uniform for a particular initial condition, there 
is no assurance that the flow can turn uniformly to a parallel flow at the exit 
for any initial condition. Recompression (i.e., deceleration of the flow) may 
occur, as indeed it does for the present initial conditions used. This recom- 
pression is avoided by arbitrarily taking the nozzle exit to be a short distance 
upstream from where the nozzle wall becomes horizontal. We thus select the exit 
to be at x = 8 cm. 

The numerical method of computation is the well-known method of characteris- 
tics. The necessary basic relations for the computation (eqs. (ill-l)) have 
already been developed in chapter III. In finite -difference form, these equa- 
tions are (with reference to fig. IV-3) 


(3 fAD ( p D - P A ) + (p w2 ) AD (e D - 9 a ) + pw2 sin 0 - A- A AD = 0 

h (iV-la) 

P %D (p D " V - (pw2) BD (0 D “ V + (l^T PW2 Sln 6 ' A]3I) = ^ 

( p w) CD (v D - v c ) + (p D - P c ) = 0 (iv-lc) 

(a D - a c ) - (I) = 0 (iV-ld) 

\ ' CD 

where single-letter subscripts refer to locations and double -letter subscripts 
signify average values, except with A, which is the linear distance between the 
two points referred to by the subscripts. Points A, B, and C are locations at 
which the values are known. Point D is the location at which the flow proper- 
ties are to be found. Here a is the degree of dissociation of oxygen. 




Figure IV- j.- Characteristics net for nonequilibrium-flow computations 



Computed array 



At the nozzle wall, equation (iV-la) is replaced by 0p = tan"" 1 S(x-p), 
where B(x) is the slope of the nozzle wall. At the nozzle center line, 
equation (iY-lb) is replaced by 0-^ = 0. 

The numerical procedure is as follows (cf. fig. IV -4) : The values of 
p, 0, w, and a are taken to be known at each point, represented by the solid 
dots, in the initial array. We shall call those points the data points. Points 
Di, D 2 , and points of type D 3 constitute the array of points to be computed. 

For the nozzle center -line calculation, B x is located halfway between data 
point Cx and the adjacent data point. The values of p, 0, w, and a at B x 
are obtained by linear interpolation. Point D ± is located at the intersection 
of the nozzle center line and the right -running Mach line from B±. For the 
nozzle wall calculation, A 2 is located halfway between data point C 2 and its 
adjacent data point, and D 2 lies on the intersection of the nozzle wall and the 
left -running Mach line from A 2 . The values of p, 0, w, and a at A 2 are again 
obtained by linear interpolation. Finally, for a typical calculation that 
involves neither -the nozzle wall nor center line, A 3 is located halfway between 
data point C 3 and the data point below it. Points of type D 3 are then taken 
to lie on the intersection of the streamline from C 3 and the left -running Mach 
line from A 3 , and B 3 at the intersection of the right -running Mach line 
through D 3 and the straight line connecting C 3 with the data point next above 
it. The values of p, 0, w, a at A 3 and B 3 are obtained by linear 
interpolation. 

For the starting of the integration process, the values of p, 0, w, a for 
points Dj_, D 2 , and type D 3 are the same as those for points C x , D 2 , and 
type C 3 , respectively. The numerical integration for points of type D 3 
involves solving equations ( IV-la) and (iV-lb) for pp 3 and 0p 3 , then equa- 
tions (iV-lc) and (iV-ld) for w-n and a p . The new values for pp , 0p , v Doj 

O O o vl5 

and ap are used in recomputing the coefficients in the compatibility relations 
and in refining the locations of the points D 3 and B 3 . The values of p, w, 0, 
a for B 3 are then reinterpolated, and the values for p, w, 0, a at B 3 are 
recomputed. The computation process is then repeated. The iteration is con- 
tinued until the velocity reaches a preassigned accuracy. For errors of 
10 ” 4 percent (between successive iterations) or less in the velocity, it usually 
takes four iterations to compute one point. The computations for T)x and D 2 are 
essentially the same except that one less computation is needed because the value 
of 0 is known at D x and D 2 . After the entire new array consisting of 
points Di, D 2 , and type D 3 are computed, the computed array is used as a new 
initial array. 

There is very little Improvement in the accuracy of the computed results 
once the number of points in an array exceeds 18. The grid sizes, which are 
dependent on the number of data points in each array, are selected to be as small 
as necessary to give the desired accuracy for the flow quantities at the nozzle 
exit. The selected number of data points in each array in the actual computation 
is roughly 30. The streamwise mesh size for a 30 -point array is approximately 
0.03 cm at the start of the nonequilibrium calculation. This mesh size is about 
one-half of the corresponding characteristic reaction length T r w. The stream- 
wise mesh size becomes larger as the flow proceeds downstream as a result of the 
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widening of the nozzle and the increasing of the Mach number. The streamwise 
mesh size at the exit is approximately 0.1 cm. The mesh size for all four cases 
computed (equilibrium and nonequilihrium) is about the same. 

For frozen flow, equations (iV-l) are immediately applicable; the terms 
involving w in equations (lV-la,b, and d) are simply zero. In the equilibrium 
limit, the term containing go in equation (ll-4a) combines with the leading term 
to yield a quantity M e 2 dp/ds. For the equilibrium flow, the compatibility 
equations become simply 

P SAD (p D ~ P A } + (p w 2 W 0 D - 0 A) + (Vr 1 pW2 Sln °) = 0 (TV -2a.) 

\ ii e x /AD 


^ e BD^D ^B^ ~ + C M r pw2 S ^ U ^ (lV-2b) 

\ ^e- 1 / BD 

and are essentially the same as those for frozen flow except that the Mach number 
and the corresponding characteristic directions are based on the equilibrium 
speed of sound. 

The numerical computation described above was performed on the IBM 7^90 
computer of the Ames Research Center, NASA. The machine program is coded in 
FORTRAN language. Fourteen subroutines are used to facilitate both the machine 
program verification and any future modification such as using different gas 
models (as long as only one nonequilibrium process is allowed). 

The numerical examples computed for the present study are based on a 
reservoir condition of T = 9000° K and p = 1000 atmospheres. Four cases were 
studied; equilibrium flow and three cases of nonequilibrium flow with the 
effective characteristic reaction time T r e ff taken as 0.1, 1, and 10 times 

the value of r r given by equation (A-9). The case of equilibrium flow has, of 
course, an effective characteristic reaction time of zero. 

To obtain the initial conditions, we assume that the gas is in equilibrium 
up to the location at which the Mach number M^* is 1.1 and that ahead of this 
location the gas behaves like real air. This last assumption gives a starting 
condition (i.e., values of p, p, and M^) corresponding exactly to real air with 
the reservoir conditions stated above. For the present computation these start- 
ing values of p and p are obtained from the charts prepared by Yoshikawa 
(ref. 27 ) , and the corresponding value of a = a e is computed using equa- 
tion ( A-b ) . The initial conditions are thus more realistic than those that 
would have been obtained using the simplified air model throughout, since at 
temperatures higher than those at the throat the nitrogen is dissociated. The 
assumption that the flow upstream of the location of Mf = 1.1 is in equilibrium 
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is justified, since the subsequent computations show that the flow has no 
tendency to deviate from the equilibrium curve (fig. IV-5) until downstream of 
this point. 
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Figure IV-5-" Degree of dissociation along the nozzle center line in the 
transition region for various degrees of nonequilibrium and various 
assumed nonequilibrium starting positions. 

The flow for the equilibrium condition is first computed. For a nonequilib- 
rium calculation, the value of the effective characteristic reaction length 
Tr eff w from the equilibrium flow along the nozzle is then observed. Once 

T r e ff w becomes moderately long (in the present computations about twice the 

streamwise mesh size), the flow is assumed to be in a nonequilibrium condition, 
and the computation for nonequilibrium is used for the flow downstream of this 
point. To assure that the selected effective starting point for the nonequilib- 
rium calculation is a correct one, at least one more case using a starting point 
upstream or downstream of it must be computed. The degree of dissociation along 
the nozzle center line in the transition region is presented in figure IV-5* As 
seen from this figure, for each effective T r , there is a certain location ahead 
of which the nonequilibrium effect is so small that any point upstream of it can 
be used as an effective starting point without serious error in the transition 
region. The starting point cannot be too far upstream, however, since too high 
a reaction rate would make the numerical computation unstable. 

The characteristic reaction length T r w for the equilibrium flow and the 
actual nonequilibrium flow (effective r T = 1.0T r ) is presented in figure IV -6. 

As seen from this, the characteristic length T r w increases rapidly in the 
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Figure IV -6.- Characteristic reaction length in the transition region of the 

nozzle . 



region immediately after the throat, causing the flow to deviate from equilibrium. 
Later T r w increases to values even larger than the physical length of the 
nozzle, and the flow becomes nearly frozen shortly past the conical section. 

The variation of the flow quantities along the center line are presented in 
figures IV -7 through IV -11. These variations are qualitatively the same as those 
from calculations of quasi -one -dimensional flow (refs. 5, 8; and 9)- Thus, 
quasi -one -dimensional -flow computations can give an assessment of the effects of 
nonequilibrium, at least qualitatively. 

Note that while all flow quantities along the center line deviate from the 
equilibrium-flow values, the extent of the deviation depending on the effective 
reaction rate, the deviation of the velocity (fig. IV- -8) is relatively slight. 

This phenomenon may be useful for some approximations. The temperature 

(fig. IV -10), on the other hand, deviates considerably from the equilibrium -flow 

value, causing the Mach number (fig. IV -11) to deviate correspondingly. 



Axial distance, cm 


Figure IV-7-- Degree of dissociation along the nozzle center line for various 

degrees of nonequilibrium. 
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Figure IV-10.- Temperature along the nozzle center line for various degrees 

of nonequilihrium. 
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The gas next to the nozzle wall travels a longer distance than that on the 
nozzle center line. Furthermore the flow on the center line expands more rapidly 
and hence freezes sooner. It may he conjectured, therefore, that the flow 
remains closer to equilibrium at the wall than at the center. If the flow is 
significantly out of equilibrium, then, the value of a at the wall should he 
lower than that at the center. If the flow is near equilibrium, of course, the 
value of a is lower at the nozzle center line than that at the wall simply 
because the flow along the center line expands more. This is shown to be the 
case in figure IV - 7 * 

The variations of the flow properties across the nozzle -exit plane are 
presented in figures IV -12 through IV -l6. These figures show that, in general 
(except a), the nonuniformities become more pronounced as the effective value 
of T r increases, that is, as the flow tends toward the frozen condition. 

The velocity, fairly uniform when the flow is in equilibrium (fig. IV- 13 )> 
stays fairly uniform when the flow is out of equilibrium. Note that, although 
a varies across the exit plane (fig. IV- 12 ), the corresponding values of 
Z = 1 + (a/4. 764), which dictates the relative importance of dissociation in the 
thermal equation of state, are actually fairly constant. Pressure (fig. IV- 14 ) , 
on the other hand, has a deviation from the center to the nozzle wall some 
30 percent (when effective T r = 1.0r r ) more than the equilibrium case. 

We may therefore conjecture that, since we cannot in the design of a nozzle 
by the inverse method have all of the exit quantities exactly uniform when the 
flow is out of equilibrium, we will attain the over -all best uniformity by taking 
the pressure uniform, the streamlines parallel (0=0), and the velocity and a 
unspecified. We may conjecture, also, that for a nozzle designed to give a 
uniform flow at the exit for equilibrium conditions, the pressure and Mach number 
may vary significantly from the nozzle center to the wall if the flow turns out 
to be nonequilibrium. 



Figure IV- 12 .- Variation of degree of Figure IV-I3.- Velocity distribution 
dissociation across the nozzle- across the nozzle -exit plane, 

exit plane. 
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Figure IV -14. - Pressure distribution 
across the nozzle -exit plane. 


Figure IV-15.- Temperature distribution 
across the nozzle -exit plane. 



Figure IV-16.- Frozen Mach number 
distribution across the nozzle - 
exit plane. 


In the present case, the correct 
nonequilibrium flow ( T r e ff = T r) is 

nearly frozen shortly after it passes 
the conical section, which here is a 
relatively small portion of the nozzle. 
If a longer conical section were pre- 
scribed, the flow would be nearly frozen 
at the end of the cone . Then one would 
compute the conical part by quasi-one- 
dimensional flow, which is exact in this 
region, and determine the flow in the 
rest of the nozzle by a frozen, axisym- 
metric computation for a thermally per- 
fect gas (for which a and hence z 
are constants). Such calculations are 
readily made since the methods for both 
the nonequilibrium quasi -one -dimensional 
flow and the thermally perfect axisym- 
metric flow have already been developed. 
Furthermore, if the conical section is 
large enough that effective freezing 
starts before the end of the conical 
section (so that the region in which the 


flow is not quasi -one -dimensional is essentially frozen), complete uniformity at 
the nozzle exit can be prescribed. This would provide one method of designing an 
ideal nozzle for nonequilibrium flow. 


Lastly, the useful, approximate, sudden -freeze analysis for nonequilibrium 
flow (ref. 6), which was developed for quasi -one -dimensional flow, can be 
adopted for axisymmetric flow. This can be done by switching from a computation 



of equilibrium flow to one of frozen flow when appropriate conditions are reached 
on each streamline. If freezing occurs after the flow passes the conical portion 
of the nozzle, the center line will freeze first because of more rapid expansion 
there. Since along each cross section the flow at the center tends to freeze 
first, the computation along curve AD (in fig. IV - 3 ) will reach the condition 
for sudden freezing sooner than that along CD. This should cause no difficulty, 
for one can use an equilibrium -flow computation along BD and CD together with 
a frozen -flow computation along AD without incompatability . 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Nov. 15, 19^2 


APPENDIX A 


THERMODYNAMIC AND CHEMICAL -KINETIC PROPERTIES OF A SIMPLIFIED 

AIR MODEL 


In this section we set down the thermal and caloric equations of state, the 
rate equation, and the expressions for the auxiliary properties, such as the 
speeds of sound, that appear in the coefficients of the flow equations. We will 
consider only the simplified air model described briefly in chapter IV. In this 
gas model the air is assumed to be a mixture of oxygen and nitrogen with only 
the oxygen dissociated. The translation, rotation, and vibration of each con- 
stituent are taken to be in equilibrium (i.e., they have a Maxwell -Boltzmann 
distribution), and the electronic excitation is assumed negligible. The deriva- 
tion of the properties of this gas is straightforward. Only minimum details, 
enough to indicate the assumptions involved, will be given since derivations for 
more complicated gas models are available (e.g., refs. 22, 28, and 29 ) • 

The number density (number of particles per unit mass of gas mixture) for 
the constituents in the mixture are 


n n = 

iHqCI + a) 


2 2h]q( 1 + a) 

a/b 

nRa ~ 2^(1 + a) 

where a is the ratio of the mass of oxygen in atomic form to the total mass of 
atomic and molecular oxygen, a is the ratio of the mass of nitrogen to the mass 
of the oxygen in the mixture, b is the ratio of the atomic weight of nitrogen 
to that of oxygen, and ihq is the mass of an oxygen atom. The actual values 
used for the constants are listed in table I at the end of appendix A. 

The assumption that the translation of the molecules has a Maxwell -Boltzmann 
distribution permits us to write the thermal equation of state as 

P = ^(“o + n 0 2 + n N 2 ) 
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where k is the Boltzmann constant. Substituting the relation for nQ, n 0 z > 
n^ 2 , and letting R = [l 4 - (a/b) ]k/[ 2 iiiq (1 + a)], and Z = 1 + a[l + (a/b)] -1 , 
we have 


E = ZRT 
P 


(Al) 


The specific internal energy e, with the zero energy of the mixture taken 
to he that of the atomic species at absolute zero temperature, is the sum of the 
translational, rotational, and vibrational energies of the species (since 
electronic excitation is neglected) . We thus write 


e ” S n i (^translation^ + ^rotation^ + ^vibration^ “ ^dissociation^) 


where the subscript i takes the values 1, 2, 3 corresponding to the species 
0, 0 2 , ^ 2 ; respectively. Substituting the values of n^ and denoting 
EdissociatlonjA "by ©d^ we have, for the enthalpy h = e + p/p, with 

^translation-^ = (3/2)kT and (assuming rotation fully excited) E ro tation-j_ = kT 


h = 1 1 + 1 ; i 2 + i z ) RT + R is 


&■ \ ( , 3 iv A nm ■ d J 1 


(l - oo)E v i *| 3 + ( —) 




4 - 




where E v i*b^ is given by the relation 


(A2) 


k©v 1 exp(-0v 1 /T) 
Evib± ~ 1 - exp( -© v - i /T) 


Here 0y^ is the characteristic vibrational temperature. 


The rate of change of the degree of dissociation a can be obtained by 
summing the forward and reverse reaction rates of oxygen in the presence of 
atomic oxygen and -molecular oxygen and nitrogen. Here the rate equation differs 
from equation (7.34) of reference 22 only by the additional term due to the 
presence of nitrogen as a catalytic species. The rate equation thus takes the 
form 


dec 

dt 


P 

Mq( 1 + a) 



+ 


1 


- a 


2 



+ 


1 a 

2 b 



1 


2pa 2 

’ M 0 K c (1 + a)_ 


(A3) 
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where Mq is the atomic weight of oxygen, the subscripts for k^ denote the 
particular catalytic species, and K c is the equilibrium constant. The quantity 
inside the last bracket is, from the law of mass action, zero when a has its 
equilibrium value <x e . Thus we have 


dp 


p d (t) 


exp - 


CG c 


®Do x 
u 2 ) 

T J 


(A4) 


where 


p d (t) 


(1 + a)Mo 
2 


K c exp 



(A5) 


The value of K c can be expressed by means of the partition function. In 
terms of temperature, we have 


Kc = 


2mQ /'jtmok' 

Mo" 


3/2 


e rQ Jr ( 1 


- 0 . 


exp 


v 0 ; 


(p-f-)o eJ 


(p.f . ) 


exp - 


0. 


2 el 


0 . 


D, 


0 c 


T 


(a6) 


where h is the Planck constant and 0 r is the characteristic rotational 

temperature for oxygen. The ratio of the electronic partition functions is 
assumed to be constant, consistent with the earlier assumption of neglecting 
the contribution of electronic excitation to the internal energy. The value 
of this ratio is taken to be 2k. k2. 


We use the same rate constants adopted by Hall (ref. 10). They are 

/ ©T 


k f = 3 - 56 x 10 ' 

u 2 


21 T -3/2 exp I _ 


3 D, 


0; 


T 


7 ' 


cm' 


n 3 /mole sec 


(A7a) 


k f, 


k f 0 


0 _ _ 35 T 
0 


D 0 5 


(A71) 


k f . 


7, 


% = i 
o ' 3 


(A7c) 


where ©p lias the value 59*380° K. The total forward reaction rate constant. 


that is, the quantity inside the first "bracket of equation (A3), is thus 


kf = 



(A8) 
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The reciprocal of the entire quantity multiplying the last bracket of equa- 
tion (A3) has the dimensions of time and characterizes the time it takes for the 
chemical reaction to reach equilibrium. This quantity, the characteristic 
reaction time, will be denoted as T r , that is 

T r = (1 + a)M 0 /pk f (A9) 


There remain now only the frozen speed of sound and the 

speed of sound a g . These are given by the identities (ref. 22) 


equilibrium of 


a f 


2 _ 



(A10) 


h p + hoCLep 

P - ( hp + ■h<x ct 'ep ) 


(All) 


where the subscripts p, p, a denote partial differentiations, for example, 
h p = (dh/dp) p a , a Sp = (Sa e /Sp) p . Thus they are simply related to the deriva- 
tives of h and a e . If we denote the quantity in the right-hand side of equa- 
tion (A2) by f (a, T), then hp is simply 


= Ah\ = /tft /dT\ = T 

P " W p ,a ~ W a W p ,a P Wo 


whereas hp and are 


h 


P 


T 

P 



ha = 



In a similar manner, if we denote log[p d _(T)exp( -©p^/T) ] by f 2 (T), then 
a e = (da e /dp) p and a 6p = (3a e /3p) p can be obtained by differentiating 
equation (A4) (with the aid of eq. (Al)) as 


a= 


f 2 (dT/dp) p>a - (l/p) 


[(2 - a e )/a e (l - a e )] - f £ (ST/aa) p ^ p 
2 - a e f 2 


= -(1 + f 2 T) p 


a e 


(l - a e ) [1 + (a/b)] + a e . 


-l 
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and 


-l 


ae P 2 \ ia e (l - cc e ) [1 + (a/b)] + a e J^) 

The values of a e and a e for the simplified air model are presented in 
figures A-l and A-2 for temperatures up to 8000° K and for various pressure 
levels. The variations of a f /a e and T r a e are illustrated in figures A-3 
and A -4. The maximum value of ap/a e for pure oxygen is roughly twice that 
for the air model. Pure oxygen, therefore, would have greater nonequilibrium 
effects than air, this being due to the lack of nitrogen acting as a diluent. 

The characteristic chemical reaction length for pure oxygen, on the other hand, 
is several times smaller than that for air because the nitrogen molecules are 
not as efficient a catalytic body for reaction as are oxygen molecules and atoms. 



Figure A-l.- Equilibrium degree of dissociation versus temperature for 
pressure = 0.01, 1.0, and 100 atmospheres (simplified air model). 



40 



Normalized equilibrium sonic speed, o e /a et 



0 I OOO 2000 3000 4000 5000 6000 7000 

Temperature , °K 


Figure A-2.- Equilibrium speed of sound versus temperature for pressure 
1.0, and 100 atmospheres (a e „ = a~ at standard condition = 3-320X10 4 c 




Figure A-3*- Ratio of frozen to equilibrium sonic speed for the simplified 

air model and for pure oxygen. 
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Figure A-k.- Normalized characteristic reaction length for the simplified 

air model and for pure oxygen. 



TABLE I . - PHYSICAL CONSTANTS USED FOR THE GAS MODEL 


a 

3.3 

l 

L 

0.8755 

h 

6 . 62 5 2 x 10 ~ S7 erg sec 

k 

1.38042X10" 16 erg/°K 

mo 

2 . 66 x 10 " 23 

Mo 

l6 gm/mole 

V 

59,380° K 

% 2 

113,260° K 

0 r 

r 0 2 

2.07° K 

0 V 

v 0 2 

2,230° K 

0 lr 

v n 2 

3,340° K 
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APPENDIX B 


a 


ELjr> 

b 

e 

f 

h 

k 

*f 

K c 

“0 

M 

Mq 


n,* 


P 

q. 

*i 

r 

s 

S 

t 

T 


PRINCIPAL SYMBOLS 

mass of oxygen/mass of nitrogen in gas model described in 
appendix A 

equilibrium sonic speed 
frozen sonic speed 

ratio of atomic mass of nitrogen to atomic mass of oxygen 

specific internal energy 

any flow quantity or function 

specific enthalpy; also Planck constant 

Boltzmann constant 

forward reaction rate constant 

equilibrium constant 

mass of an oxygen atom 

Mach number (frozen or equilibrium) 

atomic mass of oxygen 

particle density of species i 

pressure 

nonequilibrium parameter 

nonequilibrium parameter associated with the ith process 

radial distance from the axis of symmetry 

streamwise distance 

specific entropy 

time 

temperature 


^5 


V 


velocity 


Z compressibility factor; for the simplified air model, 

1 + a (i + I)" 

a degree of dissociation of oxygen = mass ratio of atomic oxygen to 

total amount of oxygen 

(B v/M 2 - 1 

5 slope of the boundary 

jj. Mach angle = sin 

Q streamline angle 

©D* ©r' ©v characteristic temperatures for dissociation, rotation, and 
vibrat ion , re spec t ively 

p density of the fluid 

T r characteristic reaction time (see eq. (A9)) 

w rate of change of q 
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